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The quantum pyrochlore antiferromagnet is studied by perturbative expansions and exact diagonal- 
ization of small clusters. We find that the ground state is a spin-liquid state: The spin-spin correlation 
functions decay exponentially with distance and the correlation length never exceeds the interatomic 
distance. The calculated magnetic neutron diffraction cross section is in very good agreement with 
experiments performed on Y(Sc)Mn2. The low energy excitations are singlet-singlet ones, with a 
finite spin gap. 

PACS number(s): 75.10.Jm, 75.50.Ee, 75.40.-s 



Since Anderson proposed the resonant valence bond 
(RVB) wave function for the triangular latticed, there 
has been a lot of attention focused on frustrated lat- 
tices, because they might have a spin-liquid-like ground 
state in two- or three-dimensional lattices. The frus- 
trated systems can be classified into two different cate- 
gories: structurally disordered systems and periodic lat- 
tices. Among the last one the ground state of the S = h 
quantum Heisenberg antiferromagnet on the kagome and 
pyrochlore lattices is expected to be a quantum spin liq- 
uid (QSL). The common point of these two lattices is 
the high degree of frustration as they both belong to the 
class of the "fully frustrated lattices." The classical mean 
field description indicates a pathological spectrum with 
an infinite number of zero energy modes which prevents 
any magnetic phase transition, and produces an exten- 
sive zero temperature entropyaul. The classical critical 
properties are nonuniversal for both systemsO. 

The pyrochlore lattice consists of a 3D arrangement 
of corner sharing tetrahedra (Fig. 1). All compounds 
which crystallize in the pyrochlore structure exhibit un- 
usual magnetic r-properties: Two of them, FeF3 and 
NH4Fe 2+ Fe 3+ FgH, are known to have a noncollinear long 
range ordered magnetic structure at low temperature; the 
other compounds do not undergo any phase transition, 
but many of them behave as spin glasses, although there 
is no structural disorder at all. It is remarkable that frus- 
tration in a periodic lattice may give rise to ageing and 
irreversibility so that a conpound such as Y2M02O7 has 
been considered as a "topological spin glass"0O. 

The problem of ordering, in the pyrochlore lattice 
was initiated by Andersonla who predicted that only 
long range interactions are able to stabilize a Neel-like 
ground state. More recently, mean field studiesu have 
confirmed thpgfi-oredictions; from classical Monte Carlo 
calculationsalLO'tLII it was concluded that this system does 
not order down to zero temperature, but any constraint 
will induce magnetic ordering|I3. 
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FIG. 1. Description of the pyrochlore lattice as an fee 
lattice of tetrahedra. Solid lines connect sites interacting with 
exchange J, dashed lines with exchange J'. 

The only attempt to describe the quantum S = h 
Heisenberg antiferromagneton the pyrochlore lattice has 
been done by Harris et alS3 who studied the stability of 
a dimer-type order parameter and showed that quantum 
fluctuations play a crucial role. In this Letter, we show 
that the ground state exhibits a QSL behavior: By apply- 
ing a perturbative approach to the density operator we 
show that spin-spin correlations decay exponentially with 
distance at all temperatures with a correlation length 
that never exceeds the interatomic distance. Fluctua- 
tions select collinear modes, but the amplitude of these 
modes is extremely small. Exact diagonalization of small 
clusters shows that the spectrum of the pyrochlore lattice 
looks like the kagome one, with a singlet-triplet spin gap 
and no gap for the singlet-singlet excitations. We also 
use our developpement to calculate the neutron magnetic 
cross section. The results are in very good agreement 
with previous experimental results on Y(Sc)Mn2t-3. 
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Model. -The Hamiltonian of the quantum Heisenberg 
model on the pyrochlore lattice is 

H=-jJ2 S i-Sj, (0.1) 

(i,o) 

where the summation is taken over the nearest neigh- 
bors (NN) sites. J is the negative exchange coupling. 
We describe the pyrochlore lattice as a fee Bravais lattice 
with a tetrahedral unit cell. Each cell is exactly diago- 
nalized and coupling between cells is taken into account 
perturbatively: We call J' the exchange interaction be- 
tween NN sites in different tetrahedra (Fig. 1) and we 
make an expansion in powers of A = J'/J. The Hamil- 
tonian is rewritten as H — Hq{J) + Hi(J'), where Hq 
describes the NN interactions within a tetrahedron and 
Hi the NN intertetrahedra interactions. By using this 
trick we have transformed the initial S = \ pyrochlore 
lattice into a fee lattice with 16 states per "site" (four 
spins 5). All of the thermodynamical quantities can be 
obtained from the density operator, which we calculate 
in perturbation with respect to A: 

P = e-< 3H = p Q + ... + Pn + e(X n ), (0.2) 

where 

Pn = (-1)" / / / d^.A^e^-^ 
Jo Jo Jo 

x H 1 ...e~^-^ (in)Ha H 1 e~ (i - H ° . (0.3) 

In order to evaluate the different terms of the devel- 
opment, we have derived a diagrammatic method which 
allows a systematic expansion. This method will reported 
elsewhere. All of the following quantities have been eval- 
uated analytically to the second order in A by imple- 
menting a formal program in Mathematica on a Silicon 
Graphics. For an indication, the evaluation of the spin- 
spin correlation functions from the first to the sixteenth 
neighbors took l| months of CPU time. The method 
was tested on the spin- g chain for which the unit cell was 
two NN sites and the development was made to fourth 
order in A. Using the analytical density operator and 
fixing A = 1, we obtained quite satisfactory-jresults, the 
difference with those of Bonner and FishercS being less 
then 2 % at low temperature. Furthermore, our method 
is more efficient when the characteristic length of the 
ground state is short, as expected in a spin liquid. Thus, 
the deduced characteristic length will be a control pa- 
rameter of our calculations. Moreover, contrary to usual 
high temperature expansion, our method yields to a non- 
divergent expansion when T — -» 0. 

Spin-spin correlations. -Let us consider a reference site 
So on the lattice. We define (So-Sd) — Cd as the corre- 
lation function between this site and a site at a distance 
d in the lattice. This function is easily evaluated as 

C d = ^S P [S .S d exp(-f3H)}, (0.4) 
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FIG. 2. (a) Absolute value of the spin-spin correlatin func- 
tions as a function of distance at T — 2\J\. (b) Correlation 
length deduced from the analysis of the (a)-type graphs at all 
temperatures. 

where Z is the partition function. The development was 
made up to second order in A, which allows us to calcu- 
late Cd up to the 16th neighbors. At any temperature 
between T = 10\J\ and T = |J|/100, we find that these 
correlations oscillate in sign and decay exponentially with 
the distance [Fig. 2(a)]. From these curves, we can ex- 
tract a correlation length £ defined as \Cd\ oc exp{—d/£). 
This length is temperature dependent and never exceeds 
one interatomic distance down to zero temperature [Fig. 
2(b)]; the broad feature at T w J shows the limitations 
of the method in this energy scale. 

Assuming that £ could be underevaluated for T < J 
by our method, we extrapolated £(T) from the T > J 
regime to the low T regime. Whatever the law we 
used (algebraic, exponential, stretched exponential), we 
found a saturation of the correlation length around the 
interatomic distance. Thus, our calculations are self- 
consistently controlled as the value of £ at T = OK is 
much smaller than the real spatial extension of our de- 
veloppement (five interatomic distances). 

The total spin can be obtained from the spin-spin cor- 
relation functions as 

N 

(($2 s i) 2 ) = NS ( S + 1) + 6NC i + 12NC 2 + - (0-5) 

i=l 

Since we have shown that these correlation functions are 
rapidly decreasing with distance, we can retain only a 
few terms in this expansion. Cutting the summation to 
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the second neighbors, we verify that the ground state is 
a singlet: 



1 N 

-<£>) 2 }(T = 0)-0.02 



(0.6) 



This result shows that our method is satisfactory even at 
a very low temperature, essentially because the correla- 
tions inside each unit cell are calculated exactly, and in 
a QSL, only short range correlations are important. 

Static structure factor.-For classical spins, pthe spec- 
trum of the pyrochlore lattice is very peculiarO: It pos- 
sesses two branches of zero energy modes. In order to 
study the effect of quantum and thermal fluctuations on 
these modes, we evaluated the structure factor, 



s m ' n ( q ) = J2a 



(0.7) 



where m and n are the indices of a site in a tetrahe- 
dral unit cell [(m,n) <E {1,2,3,4} 2 ]. Ca is defined in Eq. 
(0.4), q is a vector of the first Brillouin one (BZ), and 
R™' n is the vector of length d that links the sites of type 
m and n. For each q this structure factor is a 4 x 4 matrix 
whose eigenvalues give the fluctuation modes of the sys- 
tem, the lowest energy mode corresponding to the largest 
eigenvalue ujm(q)- To the first order in A, remains 
nondispersive over the entire BZ, and degeneracy is not 
lifted. To the second order, a maximum appears on the 
axis A of the BZ (Fig. 3). This maximum corresponds 
to a collinear phase where the total spin vanishes on each 
tetrahedron (see the inset of Fig. 3), and the phase be- 
tween two neighbouring tetrahedra is equal to ir. This 
confirms the results obtained for classical spins by pre- 
vious authorsliinij. We note that the degeneracy is very 
weakly lifted (1/10 6 of the width of the spectrum), but 
this is not a numerical artifact since, in our method, the 
precision can be as small as we want. A-similar behavior 
was also observed in the kagome latticetZl. 

Neutron cross section.-The neutron magnetic cross 
section can be expressed in terms of the correlation func- 
tions as 



d 2 a 
dflduj 



{1,2,3,4} 2 



J- e-^ T ™- T ^U m , n (Q,u;), (0.8) 



where 



u mtn (Q,u>)= y\ 



-q.(Ri-R j ) p -q.(T m -T n ) 



(0.9) 



x / (S i!m (0)S jtn (t))e- iut dt. (0.10) 



T m and i?j are the translations that define the position 
of a site of type m in the unit cell i of the space group 
Fdim] Q = k + q, where k is a vector of the reciprocal 
lattice and q belongs to the first BZ. From the static cor- 
relation functions calculated above, we obtain the total 
magnetic cross section dcr/dtt. We report the results as a 
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FIG. 3. The two largest eigenvalues of the structure factor 
along the A axis. The maximum corresponds to the collinear 
phase shown in the inset. 



contour plot in the reciprocal ([00/i],[/i/i0]) plane (Fig. 4). 
The absence of a signal in the first BZ indicates that the 
ground state is a singlet. The intensity is maximum at 
about Qt = [200] ± [§§0] or Q' x = [002] ± [§§0] (Fig. 4). 
These maxima are associated with long range correlations 
describing a structure where consecutive tetrahedra are 
in phase. This is different from the result obtained above 
from the structure factor. In fact, we have also calculated 
the cross section in the ([M)0],[0/i0]) plane and found an- 
other maximum at Q2 — [210] which corresponds to a 
phase 7T between tetrahedra as expected. Comparison of 
the intensities of the peaks at Qi an d Q2 favors the sec- 
ond structure, but the difference is small. So we conclude 
that there are two characteristic modes in this system, a 
7r-dephased one and an in-phase one, with a larger weight 
for the first mode. 

It is interesting to compare our results (Fig. 4) with 
the measurements made on Y(Sc)Mn2 by Ballou et al&3: 
The experimental results correspond exactly to our re- 
sult (Fig. 4). In fact, the measurements were made for a 
given energy w, but the shape of the neutron cross section 
in reciprocal space was found to be nearly independent of 
the energycJ, as if x{l^ w ) — f{l)g{u)- Thus, integrating 
the experimental results over the energy does not change 
the q dependance and we can make a direct comparison 
between our calculations and the measured map. First 
we reproduce all the main features previously obtained, 
i.e the absence of a signal in the firt BZ associated with 
the singlet ground state and the maxima at Q\ and Q[. 
Second, we find a 7r-dephased mode in the ([/i00],[0/i0]) 
plane that was also observed in the same compound. The 
existence of these two modes could explain the first or- 
der character qf the transition observed in the pure YM112 
compouncOcJ associated with the freezing of the mode 
mixture in the 7r-dephased one. Finally, the half-width of 
each peak provides information on the correlation length: 
It is found to be of one interatomic distance in experi- 
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FIG. 4. Map of the neutron cross section in the plane 
( [OOh] , [hhO] ) . The white regions indicate the maximum of 
intensity in the q space. 



merits as well as in our calculations. All these results 
indicate that a large part of the physics of Y(Sc)Mn2 is 
related to quantum fluctuations in this highly frustrated 
structure, although the magnetism of Mn is not localized 
but itinerant. 

Low energy part of the spectrum.— We performed exact 
diagonalization of small clusters, up to 12 spins i. These 
sizes are not large enough to give quantitative conclu- 
sions, but allow a first characterization of the spectrum. 
We find that the ground state is a singlet with a spin gap 
between this singlet and the lowest triplet (A « 0.7|J|). 
Inserted in this gap, there are several singlet states whose 
number is growing with the cluster size. This indicates 
that the lowest energy excitations are singlet-singlet-like. 
A similar property was also observed in-exact diagonal- 
izations performed on the kagome latticecil. It seems that 
both systems belon to the same class of QSL while the 
ID integer spin chains, for which the lowest energy ex- 
citations are spin excitations (singlet-triplet), are differ- 
ent. In our case, the relevant parameter is the topological 
frustration, while in ID Heisenberg spin liquids the low 
dimensionality plays a crucial role. 

In conclusion, we have studied the quantum Heisen- 
berg spin-i Hamiltonian on the pyrochlore lattice. It 
appears from the spin-spin correlations and the low en- 
ergy spectrum that the ground state of this system is 
a QSL similar to the kagome lattice ground state. Us- 
ing ou results we calculated the magnetic neutron cross 
section, and reproduced almost exactly the map experi- 
mentally observed in Y(Sc)Mn2. This confirms that this 
compound is a 3D QSL. 

The pyrochlore antiferromagnet is certainly the first 
example of a 3D QSL. For this class of QSL the dimen- 
sionality of the lattice seems to play a minor role, but 
we can expect that any perturbation will deeply mod- 
ify the low energy spectrum. Thus, such a disordered 



magnetic ground state should be extremely sensitive to 
chemical disorder. In £a£t, in several compounds, such 
as M-doped Y(Sc)Mn 2 El, Y2M02O7H, and Al-doped (3- 
Mncfj, disorder, even nondetectable, transforms a QSL 
into a quantum spin glass state. All of these compounds 
have in common unconventional behaviors (i.e, the sus- 
ceptibility is spin-glass-like, while the low T specific heat 
increases as T 2 and strong fluctuations of the local mag- 
netization are observed in the "frozen" state) . A quanti- 
tative understanding of disorder effects in QSL requires 
more studies, both experimental and theoretical. 

We are grateful to R. Ballou, C. Lhuillicr, and F. Mila 
for helpful discussions and comments. 
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